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A model for the infrared sector of SU{2) Yang-Mills theory, based on magnetic vortices represented by 
(closed) random surfaces, is presented. The model quantitatively describes both confinement (including the 
finite-temperature transition to a deconfined phase) and the topological susceptibility of the Yang-Mills ensem- 
ble. A first (quenched) study of the spectrum of the Dirac operator furthermore yields a behavior for the chiral 
condensate which is compatible with results obtained in lattice gauge theory. 



Diverse nonperturbative effects characterize 
strong interaction physics. Color charge is con- 
fined, chiral symmetry is spontaneously broken, 
and the axial U{1) part of the flavor symme- 
try exhibits an anomaly. Various model ex- 
planations for these phenomena have been ad- 
vanced; to name but two widely accepted ones, 
the dual superconductor mechanism of confine- 
ment, and instanton models, which describe the 
^7^(1) anomaly and spontaneous chiral symmetry 
breaking. However, no clear picture has emerged 
which comprehensively describes infrared strong 
interaction physics within one common frame- 
work. The vortex model presented here |^,^ aims 
to bridge this gap. On the basis of a simple 
effective dynamics, it simultaneously reproduces 
the confinement properties of SU{2) Yang-Mills 
theory (including the finite-temperature decon- 
finement transition), as well as the topological 
susceptibility, which encodes the ?7a(1) anomaly. 
Furthermore, a first (quenched) study of the spec- 
trum of the Dirac operator yields a behavior for 
the chiral condensate which is compatible with 
results obtained in lattice gauge theory; this in- 
dicates that the model also correctly reproduces 
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the spontaneous breaking of chiral symmetry. 

Center vortices are closed chromomagnetic 
flux lines in three-dimensional space; thus, they 
are described by closed two-dimensional world- 
surfaces in four-dimensional space-time. In the 
SU (2) case, their magnetic flux is quantized such 
that they modify any Wilson loop by a phase fac- 
tor (—1) when they pierce an area spanned by the 
loop. To arrive at a tractable vortex model, it is 
useful to compose the vortex world-surfaces out of 
plaquettes on a hypercubic lattice. The spacing 
of this lattice is a fixed physical quantity (related 
to a thickness of the vortex fluxes), and represents 
the ultraviolet cutoff inherent in any infrared ef- 
fective framework. The model vortex surfaces are 
regarded as random surfaces, and an ensemble of 
them is generated using Monte Carlo methods. 
The corresponding weight function penalizes cur- 
vature by associating an action increment c with 
every instance of two plaquettes which are part 
of a vortex surface, but which do not lie in the 
same plane, sharing a link (note that several such 
pairs of plaquettes can occur for any given link). 

Via the definition given above, Wilson loops 
(and, in complete analogy, Polyakov loop correla- 
tors) can be evaluated in the vortex ensemble, 
and string tensions extracted. For sufficiently 
small curvature coefficient c, one finds a confined 
phase (non-zero string tension) at low tempera- 
tures, and a transition to a high-tcmperature de- 
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confined phase. For c = 0.24, the SU{2) Yang- 
MiUs relation between the deconfinement temper- 
ature and the zero-temperature string tension, 
Tc / = 0.69, is reproduced. When further- 
more setting (To = (440 MeV)^ to fix the scale, 
measurement of ctoo^ yields the lattice spacing 
a = 0.39 fm. The full temperature dependence of 
the string tensions is displayed in Fig. ^ 




Figure 1. Confinement properties as a function of 
temperature, with c = 0.24. Crosses: string ten- 
sion between static quarks; circles: spatial string 
tension cr^. Whereas the former has largely been 
fitted using the freedom in the choice of c (see 
text), the latter is predicted. In the deconfined 
regime, it begins to rise with temperature; the 
value (Ts{T = 1.67 Tc) = 1.39 ctq reproduces the 
one obtained in full SU{2) Yang-Mills theory ||]. 

Note that the confined and deconfined phases 
can alternatively be characterized by certain per- 
colation properties of the vortex clusters. 

Complementarily, also the topological proper- 
ties of the Yang-Mills ensemble encode important 
nonperturbative effects. The topological charge 
Q of a vortex surface configuration is carried by 
its singular points |^,||,||, i.e. points at which the 
set of tangent vectors to the surface configuration 
spans all four space-time directions (a simple ex- 
ample are surface self-intersection points). Since 
a vortex surface carries a field strength charac- 
terized by a nonvanishing tensor component asso- 
ciated with the two space-time directions locally 
orthogonal to the surface ^ , these singular points 
are precisely the points at which the topological 



charge density e^^^r Tr F^^F\t is non- vanishing. 
Note that, in order to carry nontrivial global 
topological charge, surfaces must be non-oriented; 
lines along which the orientation flips can be as- 
sociated with Abelian magnetic monopoles 
Generic vortex surfaces in the ensemble studied 
here are indeed non-orientable. In practice, iden- 
tifying all singular points of the hypercubic lattice 
surfaces used in the present model involves resolv- 
ing ambiguities reminiscent of those contained 
in lattice Yang-Mills link configurations. The re- 
sulting topological susceptibility x — (Q^)/^) 
where V denotes the space-time volume under 
consideration, is exhibited in Fig. ||. 




Figure 2. (Fourth root of) the topological sus- 
ceptibility as a function of temperature, with 
c — 0.24. Also this result is quantitatively com- 
patible with the one in full Yang-Mills theory 0. 

A comprehensive description of the nonpertur- 
bative phenomena which determine strong inter- 
action physics must furthermore include the cou- 
pling of the vortices to quark degrees of freedom 
and the associated spontaneous breaking of chiral 
symmetry. The latter can be quantified via the 
chiral condensate ('0V'), which is related to the 
spectral density p(A) of the Dirac operator in a 
vortex background via the Casher-Banks formula 
(^''0) = 7r/9(0). One can locally associate the chro- 
momagnetic flux represented by an arbitrary vor- 
tex surface with a continuum gauge field ^ ; this 
allows to construct the Dirac operator directly in 
the continuum, preserving exact chiral symmetry. 
Globally, one must allow for the gauge field to be 
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defined on different space-time patches, in order 
to accomodate tfie non-orientability of tlie vortex 
surfaces. Wfiile the gauge fields on each patch 
can be chosen Abelian (i.e. in 3-direction in color 
space) , the transition functions in the overlap re- 
gions in general have to be non- Abelian. Since the 
spectrum of the Dirac operator here is obtained 
using the finite element method, it is natural to 
use the different finite elements as the space-time 
patches. A sample Dirac spectrum generated us- 
ing this method is displayed in Fig. ||. 
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Figure 3. Dirac spectral density p{X) in a universe 
of volume (1.17 fm)"*, with c = 0.24. 

Note the anomalous enhancement of p(A) near 
A = 0, in accordance with the divergence ex- 
pected for a Dirac operator with good chiral prop- 
erties Also the linear extrapolation of the 
bulk of the spectrum to A = 0, i.e. the spec- 
tral density obtained by truncating the A = 
divergence, quantitatively lies in the range ob- 
tained in full SU{2) Yang-Mills theory |]. Of 
course, the chiral condensate by itself does not 
directly represent a renormalization group invari- 
ant, physical quantity, but can only interpreted in 
conjunction with a definition of the current quark 
masses. These can only be reliably fixed in the 
present model after more detailed (ultimately un- 
quenched) measurements of hadronic properties. 

It is also instructive to consider a change in 
the construction of the gauge field which trivial- 
izes the topology. This can be done by replac- 
ing the monopole loops on the vortex surfaces by 
double vortex sheets spanning areas bounded by 



the loops. Then, the surfaces can be oriented, 
the global topological charge vanishes, and pre- 
sumably also many local clusters of nontrivial 
topology are trivialized (although the topologi- 
cal density as such does not necessarily vanish). 
In the Dirac spectrum, this leads to a truncation 
of the A = divergence, and instead the finite- 
volume gap near A = familiar e.g. from chiral 
random matrix theory opens, without the bulk of 
the spectrum being qualitatively changed. While 
the A = divergence thus seems intimately re- 
lated to topological properties chiral sym- 
metry breaking as such merely appears to require 
randomness of the gauge fields in a very general 
sense, without specific reference to topology. 

Several improvements of the present model are 
envisaged. For one, the treatment must be gener- 
alized from SU{2) to SU{3) color. Also the Dirac 
operator can still be developed further to explic- 
itly take into account a finite thickness of the vor- 
tices; the above construction, while embodying 
the correct topology, still deals with magnetic flux 
localized on infinitely thin world-surfaces. Even- 
tually, it is hoped that this model will become a 
useful tool for phenomenological considerations. 
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